High-order harmonic generation from diatomic molecules with large internuclear 
distance: The effect of two-center interference 



00 

o 
o 

(N 



Y. J. Cheni, J. Lm^* 

1. Institute of Applied Physics and Computational Mathematics, P.O.Box 100088, Beijing, P. R. China 
2. Craduate School, China Academy of Engineering Physics, 
P.O. Box 8009-30, Beijing, 100088, P. R. China 
(Dated: March 1, 2008) 

In the present paper, we investigate the high-order harmonic generation (HHG) from diatomic 
molecules with large internuclear distance using a strong field approximation (SFA) model. We 
find that the hump and dip structure emerges in the plateau region of the harmonic spectrum, and 
the location of this striking structure is sensitive to the laser intensity. Our model analysis reveals 
that two-center interference as well as the interference between different recombination electron 
trajectories are responsible for the unusual enhanced or suppressed harmonic yield at a certain 
order, and these interference effects are greatly influenced by the laser parameters such as intensity. 
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I. INTRODUCTION 

The high-order harmonic generation (HHG) from 
atoms and molecules has been one of the most intensely 
studied aspects of strong-field physics [l|, 0, d, 0, d, 0, [3, 
H, H, [13, [m, [13, [H, [13, [l3| , because it can be apphed as a 
coherent ultrashort radiation source in the extreme ultra- 
violet (XUV) and soft x-ray regions [l6|. Due to the multi- 
center characteristics and greater freedom, the molecules 
show more complicated structures in their HHG spec- 
trum than atoms and the dynamics of molecules in the 
external field is more difficult to investigate. 

Am ong many treatments, the semiclassical recollision 
model[l3| is an applaudable theory that provides a ap- 
propriate picture for HHG. In this model, the high har- 
monics are generated by a three-step sequence: (1) tunnel 
ionization of the highest energy of the electron, (2) accel- 
eration of the free electron in the laser field, and driving 
the electron back to the parent ion, (3) recombination of 
the electron to the state from which it originated. The 
recombination step leads to the emission of a XUV pho- 
ton whose energy is given by the sum of the electron of 
kinetic energy plus the ionization potential of the state. 
A widely used semianalytical approach to describe HHG 
is the strong field approximation (SFA) or Lewenstein 
model 18] , which is applicable to the molecules. It can be 
regarded as the quantum mechanical version of the semi- 
classical recollision model. The Lewenstein model has 
the advantage of requiring less computational effort than 
the ab initio solution of the time-dependent Schrodinger 
equation, which becomes very demanding at high laser 
intensities. Another advantage of the model is that it 
yields a physical interpretation of the underlying mecha- 
nisms and a certain degree of analytical description. 

Recently, several theoretical and experimental reports 
study two-center interference on the HHG from molecules 
of H+[ll,[l[l[l and C02[i[E[3,l- Lein et al observed 
the phenomenon and attributed it to the returning elec- 
tron wave colliding with cores of the two-center molecule, 
where an interference occurs fl'l. They found that the 



interference effects are sensitive to the molecular orien- 
tation but not to the laser parameters. Ciappina et al 
use the SFA model to study two-center interference [Tsj. 
They compared the SFA prediction of the interference- 
minimum position to the exact value obtained from the 
ah initio calculations. They observe that using the two- 
center continuum wave functions, instead of a plane wave 
to describe the continuum electron, greatly improves the 
comparison. The usage of the continuum wave functions 
can be considered an attempt to account for Coulomb 
effects on the returning electron. 

The previous discussions, however, are mainly toward 
the molecules with small internuclear distance R. De- 
tailed investigations for the interference effects on HHG 
from molecules with large R are still in lack. For the 
large internuclear distance case, strong coupling between 
the ground state and the first excited state emerges due 
to charge- resonance effect [l3|- This could complicate the 
HHG process. 

In the present paper, using a SFA model that considers 
the Coulomb effects, we investigate the HHG from a ID 
model diatomic molecule with large R (16 a.u.). The SFA 
model calculation gives the HHG spectrum qualitatively 
consistent with the numerical result and shows the pro- 
nounced hump and dip structure in the plateau region. 
Detailed analysis shows that two-center interference as 
well as the interference between different recombination 
electron trajectories are responsible for the enhanced or 
suppressed harmonic yield at a certain order. And in 
particular, the location of the unusual structure is found 
to be sensitive to the laser intensity. 

This paper is organized as follows. In Sec. II, we present 
our analytical theory. In Sec. HI, we apply our theory to 
the HHG from a ID model diatomic molecule with large 
R. We also analyze the complicated interference patterns 
here. Sec. IV is our conclusion. 
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II. ANALYTIC THEORY 

In Ref.[l2|, we have developed a SFA model for the 
HHG from symmetric diatomic molecules, emphasizing 
the influence of the charge-resonance states that are 
strongly coupled to electromagnetic fields for the case of 
large R. And we showed there that for sufficiently large 
internuclear distances, while initially the system is in the 
ground state, the contribution to the harmonic comes 
mostly from the continuum state-ground state transition. 
For very large R with D « 0, where D = (0|p|l), |0) and 
|1) are the ground state and the first excited state of the 
unperturbed system, respectively, the harmonic formula 
along the field direction can be written as 
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whereJ„ = J^^[^^)^J„^ = J^i^), A = (m - - 
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— Eq, B — {m + ~ ^2 — -^0- Eo is the ionization 
potential of the ground state, Aq — Ee/iu. E is the field 
amplitude, e is the unit vector along the field direction, 
Lu is the field frequency. Eq. (1) is applicable for the 
large internuclear distance and the weak field as discussed 
in our previous paper. In the SFA[18|, an important 
assumption is that electrons in the continuum states can 
be treated as a free particle moving in the electric field 
without considering the Coulomb potential. Accordingly, 
the continuum wave function |p) is approximated by a 
plane wave |p) — e'P '' with the energy Ep = p^/2, where 
the binding potential is completely omitted. Here, we 
extend to consider the modification from the Comloub 
potential, thus a better approximation of the continuum 
states in the following form is exploited. 



momentum p with certain m. The two delta functions 
in Eq. (3) denote the emission of harmonics with cer- 
tain m and n in the recombination process. From the 
above explicit expression, we conclude that the recombi- 
nation electrons with diverse momenta could contribute 
to a common harmonic order, since the delta functions 
explicitly depend on the integer n. This implies that, 
not only in the ionization process the electron can ab- 
sorb TO ± 1 photons, but also in the recombination pro- 
cess the electron can absorb (positive) or emit (negative) 
n photons, which induces the final emission of m ± 1 + ?i 
photons. 

The integrand of Eq. (3) can be divided into several 
parts. The part of (0|pk)(p-e)^(pk|0) explicitly depends 
on the ionization energy Eq and the internuclear distance 
R, but is independent of the field parameters. It incor- 
porates the interference effects between cores of the two- 
center molecules 0, [2l| , and mainly is determined by the 
properties of the molecules. The part of J„ Jm explicitly 
depends on the field amplitude E and the field frequency 
uj, representing the probability amplitude for an electron 
to absorb or emit n-|-TO± 1 photons. It reflects the inter- 
action between the electron and the field and is closely re- 
lated to the field parameters. Thus, those recombination 
electrons that have the momenta satisfying the condition 
of A = or B = all could contribute to a harmonic 
order, and the weight of the contribution depends on 
not only two-center interference but also the interference 
between different recombination electron trajectories [l^. 
The latter is strongly depends on the laser intensity, as a 
result, we expect that the interference structure in HHG 
spectra could rely on the field parameters. The above 
theoretical analysis is verified by our numerical simula- 
tions as shown in following. 



III. NUMERICAL RESULTS 



(2) 



where the effective momentum pk = PvP^ + 2£^o/|p|, 
and the energy Ep ~ p^/2. This incorporates some ef- 
fects of the binding potential through its dependence on 
EoU, i, d, [ill- With the rectification of Eq. (2), the 
harmonic formula of Eq. (1) can be rewritten as 
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For comparison with ID numerical simulation, the 
model is simplified to an one dimensional model and the 
integral over the momentum is evaluated by the so-called 
pole approximation[2Q| throughout this paper. By set- 
ting A — or B = 0, which denotes the energy con- 
servation in the ionization process, one can obtain the 



The Hamiltonian of the ID 



molecule studied here is H(t)=— 



model diatomic 
z 



2dx 



, ^ — x£ sm(Lut) , where Z is the effective 
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charge, R is the internuclear separation, E is the 
amplitude of the external electric field, and lo is the 
frequency of the external field. In the paper, we adopt 
the atomic units, h = e = = 1. Calculations 
have been performed for 780 nm trapezoidally shaped 
laser pulses with a total duration of 10 optical cycles 
and linear ramps of three optical cycles. Numerically, 
the above Schrodinger equation can be solved by the 
operator-splitting method|22|. 

A typical result is presented in Fig. 1. Here we plot 
the harmonic spectra of a ID model diatomic molecule 
with Z = 1, R — 16 a.u. and Eq = 0.638 a.u. at the 
field intensity J = 5.3 x 10i3W/cm2(Figs. 1(a) and (c)) 
and / = 1 X 10"W/cm2(Figs. 1(b) and (d)). Figs. 1(a) 
and (b) are the numerical results and Figs. 1(c) and (d) 
are the analytic results calculated by Eq. (3). The min- 
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imum or dip structure is at the 17th and the 15th order 
in Figs. 1(a) and (c), respectively, and that in Figs. 1(b) 
and (d) both are at the 23th order, as indicated by the 
vertical arrows. Though there are some differences in 
the accurate positions of the minimum between the an- 
alytic and the numerical results, the dip structure can 
be distinguished in all plottings of the spectra. In addi- 
tion, a hump structure around the 11th order harmonic 
in Figs. 1(a), around the 15th order harmonic in Figs. 
1(b), around the 9th order harmonic in Figs. 1(c), and 
around the 13th order harmonic in Figs. 1(d), can be 
identified, as indicated by the horizontal arrows. 
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Figure 2: Photon-emission spectra of a symmetrical diatomic 
molecule with Z = 0.6, R = 2.5 a.u. and Eo = 0.6 a.u. at the 
field intensity 7 = 1 x 10^*"W/cm^ (Figs. 1(a) and (c)) and 
/ = 1.5 X 10^''W/cm2(Figs. 1(b) and (d)). (a) and (b): The 
numerical results; (c) and (d): the analytic results calculated 
by Eq. (3). See the context for the illumination of the arrows. 



R — 12 a.u. at different laser intensity. They exhibit the 
similar phenomena as revealed in Fig. 1. 

Fig. 3 shows the relationship of the momentum p and 
the function 



Figure 1: Photon-emission spectra of a symmetrical diatomic 
molecule with Z = 1, i? = 16 a.u. and Eq = 0.638 a.u. at the 
field intensity J = 5.3 x 10^^"W/cm^(Figs. 1(a) and (c)) and 
7 = 1 X 10"W/cm2(Figs. 1(b) and (d)). (a) and (b): The 
numerical results; (c) and (d): the analytic results calculated 
by Eq. (3). See the context for the illumination of the arrows. 



In Fig. 2, we plot the harmonic spectra of a ID model 
diatomic molecule with Z = 0.6, R = 2.5 a.u. and Eq = 
0.6 a.u. at the field intensity / = 1 x lO^^W/cm^ (Figs. 
2(a) and (c)) and / = 1.5 x 10i'^W/cm2(Figs. 2(b) and 
(d)). Figs. 2(a) and (b) are the numerical results and 
Figs. 2(c) and (d) are the analytic results calculated by 
Eq. (3). The harmonic spectra in Fig. 2 show a board 
and flat suppressed region, which is nearly spreading in 
the whole plateau region, as indicated by the horizontal 
arrows. In contrast, for the large internuclear distance 
cases as indicated by Fig. 1, the harmonic spectra show 
some kind of hump or dip structure. Especially as the 
field intensity increases, the hump structure become more 
pronounced around the 15th order harmonic in Figs. 1(b) 
with / = 1 X 10^'*W/cm^, compared to that in Figs. 1(a) 
with / — 5.3 X lO^^W/cm^. Our calculations have been 
extended to the model molecules with R = 10 a.u. and 



^i(p) = (0|pk)(p-e)2(pk|0). 
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for the chosen model molecule in Fig. 1. The arrows in 
it indicate those momenta at which the minimal extrema 
of Eq. (4) appear. The relevant momentum values are 
labelled above the arrows. The neighboring peaks and 
hollows in Fig. 3, corresponding to certain maximal and 
minimal extrema of Eq. (4), are very close. 

According to the simple picture regarding the nu- 
clei as point emitters proposed by Lein et al^, for 
molecules with symmetric initial states, destructive in- 
terference (minimum in the harmonic spectrum) occurs 
when pkRcos{9) = (2n -I- l)7r,n = 0,1,2,.... Here 
Pk = |Pk| is the effective momentum of the electron, 9 is 
the angle between the molecular axis and the laser polar- 
ization direction. Constructive interference (maximum 
in the harmonic spectrum) occurs when pkRcos{9) ~ 
2mT,n = 0, 1,2,...[1 Q. In the ID case, for R = 2.5 
a.u., the first predicted minimum is at the 13th order 
harmonic, which agrees with that in Figs. 2(a). "While 
in Figs. 2(b), the minimum shifts to the 17th order har- 
monic. For the large 7? = 16 a.u., the predicted minima 
are at the 8th order with n = 2, the 16th order with 
n = 3, and the 27th order with n = 4, etc. The maxima 
are at the 12th order with n = 3 and the 21th order with 
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Figure 3: The relationship between the function 5'i(p) and 
the momentum p for the chosen model molecule in Fig. 1. 
See the context for the illumination of the arrows and the 
numbers above the arrows. 

n = A, etc. While in the numerical cases, the pronounced 
minimum is at the 17th order in Figs. 1(a) and the 23th 
order in Figs. 1(b). The pronounced maximum is at the 
15th order in Figs. 1(b). From the above observations, 
we find that the interference patterns in the HHG usu- 
ally depend on the laser intensity. Moreover, for the large 
R case, it becomes difhcult for the simple point-emitters 
model to predict the interference patterns. 

To illuminate the hump and dip structure in Figs. 1(b), 
we investigate the physical mechanism for the emission 
of a harmonic order from molecules with large R. In 
Fig. 4, using the same laser and molecule parameters 
as in Figs. 1(b), we plot the values of S'2(p) without 
considering the interference term Eq. (4) (the red curves 
with the hollow symbol) and S{p) with considering the 
interference term Eq. (4) (the black curves with the solid 
symbol) as functions of the electron momentum p for 
certain individual harmonic orders. Here 

S2{P) = J„(^^)J™(^), (5) 

and 

^(P) = (OIpXp • e)2(p|0) J„(^:P-^) J„(P-^). (6) 

UJ UJ 

And only the contributions from positive momenta are 
shown. For clarity, we use the dotted arrows to indi- 
cate those momenta, which correspond to the first several 
minimal extrema of 5i(p) in Fig. 3. The relevant mo- 
mentum values are labelled above the arrows. We use the 
horizontal arrows in Figs. 4(a) to indicate several maxi- 
mal extrema of >5'2(p) (denoted as peaki and peak2). 



Figure 4: The values of ^(p) (the red curves with the hollow 
symbol) and S'(p) (the black curves with the solid symbol) as 
functions of the electron momentum p for certain individual 
harmonic orders. Harmonic orders are as indicated. The laser 
and molecule parameters are the same as in Figs. 1(b). Only 
the contributions from positive momenta are shown. See the 
context for the illumination of the arrows and the numbers 
above the arrows. 



From Eq. (6), the contributions to single peak of HHG 
consist of two parts, i.e., the interference term >S'i(p) and 
Bessel function term S2 (p) ■ For the 2?>th order harmonic 
at which the harmonic dip appears in Figs. 1(d), with- 
out considering the interference term S'i(p), the primary 
contributions of S'2(p) come from those momenta around 
the values of p = 0.77 a.u. and p — 1.37 a.u. correspond- 
ing to two of the minimal extrema of Si (p) , as indicated 
by the red curve with the hollow symbol in Figs. 4(a). 
When considering the interference term S'i(p), the con- 
tributions from those momenta are strongly suppressed, 
as indicated by the black curve with the solid symbol in 
Figs. 4(a). In Figs. 4(a) those momenta around the value 
of p = 1.02 a.u. give positive and negative contributions 
to the amplitude of the 23th order harmonic. Because 
the positive and negative values will cancel each other so 
the sum of these contributions is small. 

For other harmonic orders, for example, the 13th or- 
der, around which the harmonic hump appears in Figs. 
1(d). The red curve with the hollow symbol in Figs. 
4(b) shows that the primary contributions of 6*2 (p) to 
the harmonic order come from those momenta around 
the value oi p = 1.02 a.u., different from p = 0.77 a.u. 
and p = 1.37 a.u. of the 23th order case. The black 
curve with the solid symbol in Figs. 4(b) also shows that 
the suppression effect from the interference term S'i(p) 
is weaker than that of the 23th order case. The above 
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observations indicate that the interference term S'i(p) is 
mainly responsible for the formation of the hump and dip 
in Figs. 1(b). 

Fig. 4 clearly shows that the formation of single order 
harmonic is closely related to the recombination electrons 
with many different momenta p. The contribution from 
the recombination electron at one certain momentum p 
is weighted by the corresponding probability amplitude 
5(p) = S'i(p)S'2(p). In addition, from Figs. 4(a) it can 
also be concluded that the interference minimum in the 
harmonic spectrum could appear at the harmonic order 
at which the primary contributions from the term 5*2 (p) 
are strongly suppressed due to the destructive interfer- 
ence from the term 5'i(p). Since the probability ampli- 
tude S{p) depends on the field intensity, the position of 
the harmonic hump or dip depends on the intensity too. 

The interference patterns in HHG for molecules with 
large R are much more complicated than that for 
molecules with small R. The adjacent maximal and mini- 
mal cxtremum points of Si (p) correspond to constructive 
and destructive interferences. They are so close as shown 
in Fig. 3, that the constructive and destructive interfer- 
ence extrema in the plottings of 5(p) could be entangled. 
This makes the interference effects on HHG very sensi- 
tive to the field intensity. In this case, the simple point- 
emitter model is not available. Our model could provide 



a good description for the HHG of diatomic molecules 
even for the case of large internuclear separation. 



IV. CONCLUSION 

In conclusion, using a SFA model that considers 
Coulomb potential modification on the continuum wave- 
functions, we have analytically and numerically investi- 
gated the HHG from diatomic molecules with large in- 
ternuclear separation. The harmonic spectra obtained 
by the SFA model agree with the numerical simulations. 
Our model calculation reveals that the two-center inter- 
ference as well as the interference between different re- 
combination electron trajectories are responsible for the 
unusual enhanced or suppressed harmonic yield at a cer- 
tain order, and these interference effects depend on the 
laser parameters such as intensity 
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